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Abstract. 

We have investigated the current-voltage characteristics (IVCs) of intrinsic Josephson 
junctions (IJJs). Recently, it is predicted that the longitudinal plasma wave can be excited by 
the parametric resonance in IJJs. Such an excitation induces a singularity called as breakpoint 
region around switch back region in the IVC. We have succeeded in the observation of the 
breakpoint region in the IVC of the mesa with 5 IJJs at 4.2 K. Furthermore, it is found that 
the temperature dependence of the breakpoint current is in agreement with the theoretical 
prediction. This suggests that the wave number of the excited plasma wave varies with 
temperature. 



1. Introduction 

Layered high T c superconductors (HTSCs) represent a stack of alternating superconducting and 
insulating layers, which form a system of IJJs. [U [2] Coupling between layers leads to the multiple 
branch structure with hysteresis in the IVC of such stack. In the system, hysteretic behavior is 
different from the case of single Josephson junction. For example, the value of return current 
depends on from which branch the system has returned to the zero voltage state. In addition, 
it has been recently predicted that the system shows a new feature called as breakpoint region 
(BPR), on the IVC. [31 HJ [5] The BPR is characterized by the breakpoint current (BPC) Ib p and 
transition current Ij. The former is a resonance point, at which the longitudinal plasma wave 
(LPW) is created in the stack with given number and with given distribution of the rotating 
and oscillating IJJs, while the latter corresponds to the jumping point from the given branch to 
the another one (see Fig. Ufa)). However, to our knowledge no precise experiment to observe 
the BPR has been carried out so far. 

Furthermore, the problem to determine the value of the coupling parameter in the stack of 
the IJJs also attracts much interest, but there are still many questions concerning this topic. 
The new opportunity to make this problem clear gives the investigation of the breakpoint in the 
current-voltage characteristics and temperature dependence of the breakpoint current. 

In this paper we present the experimental observation of the BPR as resulting from the 
longitudinal plasma excitation in IJJs. 



2. Theoretical background 

Physical properties of the IJJs in HTSC might be understood in the framework of the capacitively 
coupled Josephson junction (CCJJ) model. [El OH] Inclusion of the diffusion current (DC) leads 
to the multiple branch structure close to the equidistant one. In such CCJJ+DC model, phase 
dynamics is described by the system of equations (9J [10] 
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superconducting layers (S-layers). Here, d\ is the phase of the order parameter in S-layer I, A z 
is the vector potential in the barrier, I is the bias current, (3 and a are the dissipation and 
coupling parameters, respectively. The CCJJ+DC model is different from the CCJJ model by 
the last term on the right hand side. This coupled Ohmic dissipation term might be derived 
by the microscopic theory jDj or phenomenologically by the inclusion of the diffusion current 
between layers. [TU] 

Fig. Ufa) shows the simulation result of the IVC of a stack of 8 IJJs for a = 1, f3 = 0.2 and 
nonperiodic boundary condition 7 = 0, where 7 = s/so = s/sn and s, so, sjy are the thicknesses 
of the middle, first, and last superconducting layers, respectively. The IVC demonstrates the 
equidistant branch structure and a breakpoint on their hysteresis branch. 

As was shown in Ref. [6j, the system of equations for the CCJJ model has a solution 
corresponding to the LPW propagating along the c axis. In this case, a frequency of the LPW 
at I = and (3 = is ui p (k) = uj p \/1 + 2a (1 — cos(k(d + s)), where u> p is the Josephson plasma 
frequency, k is wave vector of the LPW and d is the thickness of the insulating layer. The 
Josephson oscillations, which frequency is determined by the voltage value, excite the LPW by 
their periodical actions, so that at ujj = 2lv p the parametric resonance is realized and the LPW 
with definite wave number is created. O IH [5] Then, the wave number depends on the number 
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Figure 1. (Color online) (a) The simulated IVC of a stack of 8 IJJs for a = 1, /3 = 0.2 and 
7 = 0; (b) Profile of the charge oscillation on the first S-layer in the stack with 7 IJJ. The inset 
shows the charge oscillation on the first layer in the beginning of the BPR. 



of junctions in the stack, boundary conditions, dissipation and coupling parameters. Using 
Maxwell equations the charge density pi on the superconducting layer I can be expressed by 
the voltages Vi and VJ+i in the neighbor insulating layers. Solution of the system of dynamical 
equations for the gauge-invariant phase differences between S-layers gives us the voltages V\ in 
all junctions in the stack, and it allows us to investigate the time dependence of the charge on 
each S-layer. 

Fig. djb) shows the profile of the charge oscillation on the first S-layer in the stack with 7 
IJJs and the inset shows the charge oscillation on the first layer in the beginning of the BPR. 
Here, the time is normalized to l/aj p o and the charge is normalized to po = ^oVo/r^, where 
7*d is Debay screening length. In this figure, the time is calculated as i - T m Io/5I, where 
T m determines time interval for averaging, Iq is an initial value of the bias current for time 
dependence recording, 51 is a step in current. It demonstrates a sharp increase of the oscillation 
amplitude corresponding to the creation of the LPW with a wave vector k = 6tt/7. Therefore, 
we can say that the observation of the BPR is an indirect evidence of the LPW excitation. 

3. Experimental results 

3.1. Samples 

Bi2Sr2CaCu20 y single crystals with critical temperatures T c of ~ 85 K were grown by a 
conventional melting method. The grown single crystals were glued on glass substrates and 
then cleaved in air to obtain their fresh surfaces. Subsequently, Au thin films of 20 — 40 nm 
thickness were deposited on the cleaved surfaces to contact well electrically. Mesas with lateral 
areas of 9 — 64 /xm 2 , consisting of 5 — 20 intrinsic Josephson junctions, were fabricated by using 
electron beam lithography, photolithography and Ar ion milling. Then the number of IJJs in the 
mesa was determined by the etching time. Fig. [2] shows schematic view and optical photograph 
of the fabricated samples. Their IVC along the c-axis direction were measured using a three 
terminal method at the temperatures ranging from 4.2 K to T c . 




Figure 2. (Color online) (a) Schematic view of the fabricated mesa structure and (b) photograph 
of the sample. 



3.2. Current-voltage characteristics 

Fig. El shows the experimental IVCs for the BSCCO stack (sample B) with an area of 24 /xm 2 
at different temperatures: (a) 4.2 K, (b) 57 K, (c) 75 K and (d) 82 K. The stack exhibits a 
large hysteresis at low temperatures and the IVCs have a multiple branch structure. From the 
number of resistive branches it is found that this mesa consists of 8 IJJs. The critical current at 
4.2 K is 460 and the corresponding density is 1920 A/cm 2 . With the increase in temperature 
the hysteresis is decreased and disappears at temperatures near T c . In addition, the increase 




Figure 3. Typical current- voltage characteristics of a stack of IJJs (samples B) at different 
temperatures. 

in temperature leads to the appearance of the breakpoint, which is indicated by arrow on the 
outermost branch in the IVC. Similar behavior were also observed for other samples. 

3. 3. Temperature dependence of the breakpoint current 

As shown in Ref. [3], the i& p depends on the values of a and j3. Therefore, it would be interesting 
to see if the predicted features can be observed in the experiment. According to the resistively 
and capacitively shunted junction model the hysteresis is related to the dissipation parameter 
and hence the decrease in the hysteresis as can be seen in Fig. [3l corresponds to the increase 
in (3 under a constant a. Therefore, the predicted features might manifest themselves on the 
temperature dependence of the breakpoint current. 

Fig. Eta) shows the temperature dependence of the critical current I c , the breakpoint current 
Ibp and the return current I r for sample B. For this sample the obtained I C (T) significantly 
differs from the Ambegaokar-Baratoff dependence, but this may be due to the influence of flux 
trapping. In comparison with the strong T dependence of I c at low temperature, the I r hardly 
depends on T for T < 40 K, and it gradually increase with T for T > 40 K. Such dependence of 
I r is typical for IJJ stacks. [TTj 1 12j In addition, we note that the T dependence of Ib p follows that 
of I r . Similar results were also obtained for the samples with larger I c . An example is shown in 
Fig. HKb). It is found that I r and /&„ show essentially the same T dependence as that shown in 
Fig. Ha). 

In Fig. [5(a), we show the simulated temperature dependence of the breakpoint current for the 
stack with 8 IJJs at nonperiodic boundary condition 7 = for a = 1 and 3, which was obtained 
by using the Ambegaokar-Baratoff dependence of I c . The details of simulation procedure of 
the temperature dependence are presented in Ref. [13] It is found that the numerical h P {T) is 
qualitatively in agreement with the experimental one. Furthermore, according to Ref.[5J the Ib p 
is related to the creation of the LPW with different wave number k at different values of the 
dissipation parameter beta, as shown in Fig. 15(b). Thus the increase in 1^ suggests that the 



(a) 



500 



(b) 



2000 




40 60 
T(K) 



Figure 4. (Color online) Experimental temperature dependence of I c , I r and i& p for (a) sample 
B (N = 8) and (b) sample No.115 (iV=30). 




Figure 5. (Color online) (a)Temperature dependence of the breakpoint current for the stack 
with 8 IJJ at different values of coupling parameter (a = 1, 3). (b) The simulated (3 dependence 
of Ib p for the stacks with 8 and 16 IJJs for a = 3. In the figure, the wave number of the excited 
LPW is indicated. 



wave number k depends on the temperature. In order to clear this point, however, a further 
experiment would be required. 

3-4- Size dependence of the breakpoint current 

We have also study the size dependence of the breakpoint current. Fig. Eta) shows the I c and Ib p 
of the mesas, which were fabricated on the same crystal, at 77 K as a function of the mesa area 
S. One can see that the I c and Ib p increase lineally with S. The linear increase in I c indicates 
that these mesas have almost same critical current density at 77 K : they have a uniform <p 
without the in-plane distribution. Thus, we attribute the increase in with S to the increase 



in I c and hence can say that the breakpoint current is an intrinsic property of IJJs because the 
value of hp/Ic is independent of S as shown in Fig. E^b) but depends on T. 
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Figure 6. (a) The S dependence of I c and Ib p of the mesas fabricated on the same crystal, at 
77 K. (b) The Ibp/Ic ratio as a function of S. 



3.5. Manifestation of the breakpoint region 

It is usually difficult to observe the BPR on the IVC because the width of BPR, defined as the 
difference between J& p and Ij is very small and hence in the most case, the Ib p coincides with Ij 
due to thermal fluctuation. As mentioned above, however, the observation of the BPR would 
strongly support the excitation of the LPW in IJ J stack. We succeeded in the observation of the 
BPR by carefully measuring the IVC. The observed BPR is shown in Fig. [71 The temperature 
is 4.2 K and this mesa has an area of 25 ^m 2 and 5 IJJs. In Fig. [7_tb), we can see clear BPR on 
the outermost branch. The width of the BPR is about 0.5 //A, and it is found to be considerably 
smaller compared with the I c value of 240 /xA. 
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Figure 7. The IVC of the mesa with 5 IJJs at 4.2 K. (a) Overall IVC. (b) IVC in the low bias 
region. 



4. Conclusion 

The current- voltage characteristics of intrinsic Josephson junctions have been investigated. Clear 
breakpoint and breakpoint region have been observed on the current-voltage characteristics. The 
temperature dependence of the breakpoint current is qualitatively in agreement with theoretical 
prediction. Our results imply that the longitudinal plasma wave can be excited in the intrinsic 
Josephson junctions. 
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